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Output Feedback Tracking Control for a Class of Uncertain Systems 
subject to Unmodeled Dynamics and Delay at Input 

Quan Quan, Hai Lin, Kai-Yuan Cai 
Abstract 

Besides parametric uncertainties and disturbances, the unmodeled dynamics and time delay 
at the input are often present in practical systems, which cannot be ignored in some cases. This 
On) ■ paper aims to solve output feedback tracking control problem for a class of nonlinear uncertain 

■ systems subject to unmodeled high-frequency gains and time delay at the input. By the additive 

CN '■ 

state decomposition, the uncertain system is transformed to an uncertainty-free system, where the 
. uncertainties, disturbance and effect of unmodeled dynamics plus time delay are lumped into a new 



X 



disturbance at the output. Sequently, additive state decomposition is used to decompose the trans- 
formed system, which simplifies the tracking controller design. To demonstrate the effectiveness, 
the proposed control scheme is applied to three benchmark examples. 
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_ Additive state decomposition, tracking, input delay, unmodeled dynamics, output feedback, 
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I. Introduction 

Tracking control of an uncertain system is a challenging problem. Most of research mainly 
focuses on systems subject to parametric uncertainties and additive disturbances [Hl-llll. Also, 
some research focuses on systems subject to uncertainties at the input, such as backlash, dead 



^ ! zone or other nonlinearities [|4l-[|5l. It is well known that unmodeled dynamics and time delay 

at the input are also often present in practical systems. For example, the unmodeled dynamics 
and time delay at the input often exist in flight control systems [Ill-jHl. These uncertainties 
at the input may produce a significant degradation in the tracking performance or even cause 
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instability if not dealt with properly. In the literature, there are some academic examples to 
demonstrate that uncertainties at the input cannot be ignored in some cases. For example, in 
||9l , the authors constructed a simple example, later known as Rohrs' example, to show that 
conventional adaptive control algorithms lose their robustness in the presence of unmodeled 
dynamics. Also, some control algorithms may lose their robustness in the presence of input 
delay, see for example the repetitive control example considered in [ITOll . Therefore, it is 
important to explicitly consider unmodeled dynamics and time delay at the input in the 
controller design. 

In this paper, the output feedback tracking control problem is investigated for a class 
of single-input single-output (SISO) nonlinear systems subject to mismatching parametric 
uncertainty, mismatching additive disturbances, unmodeled high-frequency gains and time 
delay at the input. Before introducing our main idea, some accepted control methods in 
the literature to handle uncertainties are briefly reviewed. A nature way is to estimate all 
of the unknown parameters, then compensate for them. In IfTTI . the tracking problem for a 
linear system subject to unknown parameters and the unknown input delay was considered, 
where both the parameters and input delay were estimated by the proposed method. However, 
this method cannot handle unparameterized uncertainties such as unmodeled high-frequency 
gains. The second way is to design adaptive control with robustness against unmodeled 
dynamics and time delay at the input. In [|T2ll . the Rohrs' example and the two-cart example, 
which are tracking problems for uncertain linear systems subject to unmodeled dynamics and 
time delay at the input respectively, were revisited by the Ci adaptive control. In [13|, the 
authors analyzed that their proposed method is robust against time delay at the input. The 
third way is to convert a tracking problem to a stabilization problem by the idea of internal 
model principle [14], if disturbances or desired trajectories are limited to a special case. In 
ifTSl . the problem of set point output tracking of an uncertain linear system with multiple 
delays in both the state and control vectors was considered. There also exist other methods to 
handle uncertainties. However, some of them such as high-gain feedback cannot be applied 
to the considered system directly as they rely on rapid changing control signal to attenuate 
uncertainties and disturbance. After passing unmodeled high-frequency gains or time delay 
at the input, the rapid changing control signal will be distorted a lot which will affect the 
feedback and then may destabilize the system. This explains why high-gain feedback is often 
avoided in practice. 

Compared with these existing literature, the problem studied in this paper is more general 
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since not only the uncertainties at the input but also the output feedback and mismatching are 
considered. For output feedback, the state needs to be estimated which is difficult mainly due 
to the uncertainties and disturbances in the state equation. Even if parameters and disturbance 
can be estimated, it is also difficult to compensate for mismatching uncertain parameters and 
disturbance directly. To tackle these difficulties, two new mechanisms are adopted in this 
paper. First, the input is redefined to make it smooth and bounded to handle uncertainties 
at input. As a consequence, the effect of unmodeled high-frequency gains and time delay 
at the input is always bounded. And then, to handle estimate and mismatching problem, 
the input-redefinition system is transformed to an uncertainty-free system, which is proved 
to be input-output equivalent with the aid of the additive state decompositioixi fT6J. All 
mismatching uncertainties, mismatching disturbance and effect of unmodeled dynamics plus 
time delay are lumped into a new disturbance at the output. An observer is then designed 
for the transformed system to estimate the new state and the new disturbance. Next, the 
transformed system is 'additively' decomposed into two independent subsystems in charge 
of corresponding subtasks, namely the tracking (including rejection) subtask and the input- 
realization subtask. Then one can design controller for each subtask respectively, and finally 
combines them to achieve the original control task. Three benchmark examples are given to 
demonstrate the effectiveness of the proposed control scheme. 

The additive state decomposition is a decomposition scheme also proposed in our previous 
work [fTTl . where the additive state decomposition is used to transform output feedback 
tracking control for systems with measurable nonlinearities and unknown disturbances and 
then to decompose it into three simpler problems. This hence makes a challenging control 
problem tractable. In this paper, a different control problem is investigated by using additive 
state decomposition. Correspondingly, the transform and decomposition are different. The 
major contributions of this paper are: i) a tracking control scheme proposed to handle 
mismatching parametric uncertainty, mismatching additive disturbances, unmodeled high- 
frequency gains and time delay at the input; ii) a model transform proposed to lump various 
uncertainties together; iii) additive state decomposition in the controller design, especially in 
how to handle saturation term. 

This paper is organized as follows. In Section II, the problem formulation is given and 

'in this paper we have replaced the term "additive decomposition" in fT6\ with the more descriptive term "additive state 
decomposition". 
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the additive state decomposition is introduced briefly first. In Section III, input is redefined 
and the input-redefinition system is transformed to an uncertainty-free system in sense of 
input-output equivalence. Sequently, controller design is given in Section IV. In Section V, 
two-cart example is revisited by the proposed control scheme. Section VI concludes this 
paper. 

II. Problem Formulation and Additive State Decomposition 
A. Problem Formulation 

Consider a class of SISO nonlinear systems as follows: 

X = f (t, X, 0) + hu^ + d^x (0) = Xq 

y = c^x. (1) 

Here b E M" and c G M" are constant vectors, 9 (t) E M™ belongs to a given compact set 
n CR"^^ x{t) eW is the state vector, y (t) E R is the output, d{t) E M" is a bounded 
disturbance vector, and (t) G M is the control subject to an unmodeled high-frequency 
gain and a time delay as follows: 

(s) = H (s) e-^'u (s) (2) 

where H (s) is an unknown stable proper transfer function with if (0) = 1 representing the 
unmodeled high-frequency gain at the input and r G R is the input delay. It is assumed that 
only y is available from measurement. The desired trajectory r [t) G M is known a priori, 
t > 0. In the following, for convenience, we will drop the notation t except when necessary 
for clarity. 

For system ([T]), the following assumptions are made. 

Assumption 1. The function / : [0, oo) x M" x R"' ^ R" satisfies / (t, 0, 6) = 0, and 
is bounded when x is bounded on [0,oo). Moreover, for given 6 E ^l, there exist positive 
definite matrices P E R"''" and Q E R"''" such that 

PdJ {t, X, 6) + dlf {t, x,d)P< -Q, Vx G R^ (3) 

where d^f = %E R"^''. 

Definition 1 IfTSl . The £i gain of a stable proper SISO system is defined ||G'(s)||^ = 

/•oo ^ 

/ l^f (t)|dt, where g (t) is the impulse response of G (s). 
Jo 
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Assumption 2. There exists a known stable proper transfer function C (s) with C (0) = 1 
such that ||C (s) {H (s) — l)!!^^ < eh, \\sC (s)||_^^ < e^, where Eh, eM. are positive real. 

Under Assumptions 1-2, the objective here is to design a tracking controller u such that 
y ^ r with a good tracking accuracy, i.e., y — r is ultimately bounded by a small value. 

Remark 1. From Assumption 1, since / {t, x, 6) = d^f (t, x + /is, 6) x, fi E (0, 1) by the 
Taylor expansion. Consequently, the system x = f {t, x, 6) is exponentially stable by ([3]). In 
practice, many systems are stable themselves or they can be stabilized by output feedback 
control. The following three benchmark systems all satisfy Assumption 1. 

Example 1 (Rohrs' Example). Consider the Rohrs' example system as follows [|9l|: 

. ^ 2 229 , , 

V (s) = ^ u (s) . (4) 

^ ^ ^ s + 1 s2 + 30s + 229 ^ ^ ^ ^ 

The nominal system is assumed to be y (s) = (s) here. In this case, the system dH) can 

be formulated into ([T]) as 

x = -[3 + e)x + 2u^,x{0) = 1 

y = x (5) 

where the parameter ^ = — 2 is assumed unknown and H (s) = 72:p^j:229' = 0. It is easy 
to see that Assumption 1 is satisfied. Choose C (s) = 2^^- Then Assumption 2 is satisfied 
with Eh = 0.12 and = 1. 

Example 2 (Nonlinear). Consider a simple nonlinear system as follows fT9|: 

X = -x - {1 + 9)x^ + u{t -t) + d,x{0) = 1 

y = x, (6) 

where x,y,u,d E M, the parameter 6 (t) = 0.2 sin (0. It + 1), the input delay r = 0.1 and 
d (t) = 0.5 sin {0.2t) are assumed unknown. The system ^ can be formulated into © with 
f{t,x,e) = -X - {l + e)x^ and H {s) = l,r = 0.1. It is easy to verify d^f{t,x,e) = 
— 1 — 3 (1 + 6') < —1. Therefore, Assumption 1 is satisfied. Let C (s) = Then 
Assumption 2 is satisfied with eh = and Er = 1- 

Remark 2. The Rohrs' example system in Example 1 is proposed to demonstrate that 
conventional adaptive control algorithms developed at that time lose their robustness in the 
presence of unmodeled dynamics |[9l. For the tracking problem in Example 2, there exist 
robustness issues by using exact feedback linearization [|T9l . Compared with the system in 
|fT9l , the input delay is added in Q to make system worse. The two benchmark examples 
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tell us that the uncertainties either on the system parameters or at the input cannot be ignored 
in practice when design a tracking controller, even if the original systems are stable. This is 
also the initial motivation of this paper. 

B. Additive State Decomposition 

In order to make the paper self-contained, additive state decomposition lfT6]| is introduced 
briefly here. Consider the following 'original' system: 

f {t,x,x) = 0,x{0) = xo (7) 

where x E M". We first bring in a 'primary' system having the same dimension as (|7]), 
according to: 

fp (^) -^p) 0, Xp (0) Xpfl (8) 

where Xp E M". From the original system (|7]) and the primary system ([8]) we derive the 
following 'secondary' system: 

/ (t, X, x) - fp (t, Xp, Xp) =0,x (0) = xq (9) 

where Xp E M" is given by the primary system ([8]). Define a new variable Xs E M" as follows: 

Then the secondary system © can be further written as follows: 

/ (ti i-s ~l~ i-pi -^s ~l~ -^p) fp (fi ^p: -^p) 0, Xg (0) Xq Xp^Q. (11) 

From the definition (flOl) . we have 

x{t) = Xp{t) + Xg{t) ,t>Q. (12) 

Remark 3. By the additive state decomposition, the system ([7]) is decomposed into two 
subsystems with the same dimension as the original system. In this sense our decomposition 
is "additive". In addition, this decomposition is with respect to state. So, we call it "additive 
state decomposition". 

As a special case of (|7]), a class of differential dynamic systems is considered as follows: 

X = f{t,x),x{0) = Xo, 

y = h{t,x) (13) 
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where x G M" and y G W^. Two systems, denoted by the primary system and (derived) 
secondary system respectively, are defined as follows: 

Xp fp (t, Xp) , Xp (0) Xpfl 

Vp = hp (t, Xp) (14) 



and 



Xg / (ij Xp -\- Xg) fp (t, Xp) , Xs (0) Xq Xp^Q, 

Us = h(t,Xp + Xs) - hp{t,Xp) (15) 



where Xg = x — Xp and ys = y — yp. The secondary system (|T5I) is determined by the original 
system (fT3l) and the primary system (fT4l) . From the definition, we have 



X (t) = Xpit) + Xs (t) ,y (t) = ypit) + ys{t) ,t > 0. (16) 

III. Input Redefinition and Model Transformation 

Since H (s) is the unmodeled high-frequency gain and r is the input delay, the control 
signal should be smooth (low-frequency signal) so that it will maintain its original form as 
far as possible after passing H (s) e^^*. Otherwise, the control signal will be distorted a lot. 
This explains why high-gain feedback in practice is often avoided. For such a purpose, the 
input is redefined to make control signal smooth and bounded first. This makes the effect of 
H (s) e^'^^ under control, i.e., the effect will be predicted and bounded. 

A. Input Redefinition 

Redefine the input as follows: 

u{s) = C{s) [aa{v) is)] 

where f G M is the redefined control input and ctq : M — )■ [— a, a] is a saturation function 
defined as cTq (x) =sign(x) min (|x| , a). Then is written as 

u^{s) = H{s) e-"C(s) [aa (v) (s)] 

= Cis)[aaiv)is)]+as) (17) 

where ^ (s) = C (s) (H (s) e^''^ — 1) [aa {v) (s)] represents the effect of the unmodeled high- 
frequency gain and the time delay. The function ^ (s) can be further written as 

i{s) = C {s) {H (s) - 1) e-" [aa (v) (s)] + C {s) (e^^^ - l) [a^ {v) {s)] . (18) 
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From the definition of cTq, we have snp_^^^^^ \aa {x)\ < a. In this paper C denotes the 
inverse Laplace transform. By Assumption 2, ^ is bounded as follows: 

sup le it)\ < \\C is) {H is) - 1)11^^ a + \\sC {s)\\^^ sup |£-^ { (e^ - l) /s [a^ (v) {s)]}\ 

< e//a + sup / aa{v{X))dX 

t>o Jt 

< (eh + TEr) a (19) 

where ^ (t) = ((^ (s)) .The input redefinition makes ^ bounded not matter what the 
redefined control input v is. Therefore, the redefined control input v can be designed freely. 
According to input redefinition above, the controller ^ is rewritten as 

u^ = u + C (20) 
Here u (t) = C^^ (C (s) [aa (v) (s)]) can be written in the form of state equation as follows 

u = c^z + d^aa (v) (21) 

where the vectors and matrices are compatibly dimensioned depending on C (s) . Substituting 
(l20l) into the system ([T]) results in 

X = f {t, X, 9) + bu + dh, X (0) = Xq 

y = c^x (22) 



where dh = d + ^. The system (1221) with the redefined controller (1211) is called as the input- 
redefinition system here. 

B. Model Transformation 

The unknown parameter 6 and the unknown disturbances d are not appear in "matching" 
positions for the control input, i.e., 6 and d do not appear like b [u^ + 9^x + d) . Therefore, 
in a general system except for one dimensional system, the unknown uncertainties cannot be 
often compensated for directly. Even if 9 and d satisfy the "matching condition", it is also 
difficult to compensate for since the state x is unknown. To tackle this difficulty, we first 
transform the input-redefinition system (l22l) to an uncertainty-free system, which is proved 
to be input-output equivalent with the aid of the additive state decomposition as stated in 
Theorem 1. Before proving the theorem, the following lemma is needed. 
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Lemma 1. Consider the following system 

x = f{t,x + z,9)-f{t,z,e) + p (23) 
where p{t) G M" is bounded. Under Assumption 1, the solutions of (1231) satisfy 

</3(||x(to)||,t-to)+7 sup \\p{s)\\ (24) 

to<S<t 

where /3 is a class /C£ function f20', p. 144] and 7 = 2-;-^^f^\^. 

'^min V -'^^y '^min V ^ J 

Proof. By the Taylor expansion, the function f {t,x + z, 0) can be written as 

f {t,x + z,9) = f (t, z, 6) + dxf {t,x + z + px, 0) X 
where p E (0, 1) . Then the system (l23l) can be rewritten as 

i = dxf {t,x + z + px, 9) X + p. (25) 

Choose Lyapunov function V = x^Px. By Assumption 1, the derivative of V along (l25l) 
satisfies 

< -Amin (<5) Ikf + Amax (^) I^H ||p|| 

<-^A_(Q) \\x\\\ ^M>2^^^^\\p\\. 

^ ^min 

By Theorem 4.19 [|20l p. 176], we can conclude this proof. □ 
With Lemma 1 in hand, we have 

Theorem 1. Under Assumption 1, there always exists an estimate of 6, namely 6 E ^l, 
such that the system (l22l) is input-output equivalent to the following system: 

inew f(j^j ■^newi ^) ~l~ •^new (0) 

y C Xneui ~l~ dnew (26) 

Here x„e«) and dnew satisfy 

\\X-Xnew\\ < /^(Ikoll , ^ - ^o) + 7 SUp ||rfg(s)|| 

to<s<t 

||c^ne«;|| < ||c|| /?(||Xo|| ,t - to) + ||c|| 7 sup ||dg(s)|| (27) 

where /3 is a class JCC function, 7 = 2 ^ . '('pfA^'' (q) ^'^d dg = f (t, x^eu-, 6') — /(t, x„e«,, 0)+dh. 
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Proof. In the following, additive state decomposition is utilized to decompose the system 
(l22l) first. Consider the system (|22l) as the original system and choose the primary system as 
follows: 

Xp = f{t, Xp, 9) + hu, Xp (0) = 

Vp = c^Xp. (28) 

Then the secondary system is determined by the original system (|22l) and the primary system 
(|281) with the rule ([HI) that 

Xs = f (t, Xp + Xs, 9) - f{t, Xp, 9) + dh, Xs (0) = xq 

Us = c^Xs- (29) 



According to (fT6l) . we have x = Xp + Xg and y = yp + yg. Consequently, we can get an 
uncertainty-free system as follows 

ip = f(ti ^) + ^p (0) = 
y = c^Xp + ys 

where u and y are the same to those in (l22l) . Let and dnew = We can conclude 

that the system (l22l) is input-output equivalent to (l26l) . Next, we will prove that (1271) is 
satisfied. The system (l29l) can be rewritten as 

is = / (t, Xp + x^, 6') - / (t, Xp, 9) + rfg, (0) = xq 

Vs = c^Xs (30) 

where dg = f {t, Xp, 9) — f{t, Xp, 0) + d^. Then, by Lemma 1, we have 

\\X (t) - Xnew {t)\\ = \\Xs (t)\\ < (3 {\\xo\\ , t - to) + 7 SUp \\dg 

to<s<t 

\\dnew {t)\\ < \\c\\ \\Xs < || c|| /3 ( || Xq || ,t-to) + \\c\\ 7 SUp (s)|| . 

to<s<t 

□ 

For the uncertainty-free transformed system (l26l) . we design an observer to estimate Xnew 
and dnew, which is stated in Theorem 2. 

Theorem 2. Under Assumption 1, an observer is designed to estimate state Xnew and d^ew 
in (|26|) as follows 

Xjiew) f(j^j Xneui) ^) ~l~ ^new (0) 

dnew 1/ C Xfiew (31) 
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Then X^^-iij — ^new dnew — ^new • 

Proof. Subtracting (|3T]) from (l26l) results in 



•^new dxf(tj Xnew ~\~ ^new ~\~ fJ'-^newy ^^-^newi ■^new (0) 0, 



where ^ G (0, 1) and Xnew = Xnew - Xnew Then x^ew = 0. This implies that Xnew = Xnew 
Consequently, by the relation y = c^Xnew + dnew in (l26l) . we have dnew = dmw D 

Remark 4. By (I2TI) . the control signal is always bounded. Therefore, by Lemma 1, the 
state Xnew is always bounded. Consequently, by (l27l) . fi„eu; is always bounded as well. It is 
interesting to note that the new state Xnew and disturbance dmw in the transformed system 
(l26l) can be observed directly rather than asymptotically or exponentially. This will facilitate 
the analysis and design later. 

Example 3 {Rohrs' Example, Example 1 Continued). According to input redefinition above, 
the Rohrs' example system ([5]) can be rewritten as follows: 

x = - (3 + e)x + 2u+{d + 2^) 
y = x 

where sup4>o \C (01 < 0.12a, and u is generated by i = —0.52 + O.Sua (v) ,u = z. Then, 
according to (|26|) . the uncertainty-free transformed system of ([5]) is 

y Xnew ~l~ dnew (32) 

where 9 will be specified later. 

Example 4 (Nonlinear, Example 2 Continued). According to input redefinition above, the 
nonlinear system Q can be rewritten as follows: 

x = -x-{l + 9)x'^ + u + {d + ^) 
y = x 

where sup4>o \^ {t)\ < 0.1a, and u is generated by i = —0.52; + O.Scr^ (f ) , m = z. Then, 
according to (l26l) . the uncertainty-free transformed system of ^ is 

Xnew Xnew ~l~ ^^Xnew ~^ Xnew (0) 

y Xnew ~l~ dnew (33) 

where 9 will be specified later. 
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IV. Controller Design 



In this section, the transformed system (1261) is 'additively' decomposed into two indepen- 
dent subsystems in charge of corresponding subtasks. Then one can design controller for 
each subtask respectively, and finally combines them to achieve the original control task. 



A. Additive State Decomposition of Transformed System 

Currently, based on the new transformed system (l26l) . the objective is to design a tracking 
controller u such that y r with a good tracking accuracy, i.e., y — r is ultimately bounded 
by a small value. While, u is realized by (|2TI) . According to this fact, the transformed system 
(|26|) is 'additively' decomposed into two independent subsystems in charge of corresponding 
subtasks, namely the tracking (including rejection) subtask and the input-realization subtask. 
This is shown in Fig.l. 



Tracking (Rejection) Task with 
Saturation (26): j; — >• r 



r 



1 



Tracking (Rejection) Subtask 
without Saturation (34): — > r 



Input- Relization Subtask (35): 

y. ^ 



Fig. L Additive state decomposition flow 

Consider the transformed system (l26l) as the original system. According to the principle 
above, we choose the primary system as follows: 

Xp = f{t, Xp, 9) + hup, Xp (0) = 

Dp C Xp -\- dnew (34) 

Then the secondary system is determined by the original system (l26l) and the primary system 
(l34l) with the rule (fT5l) . and we can obtain that 

is = f{t, Xp + Xs, 9) - f{t, Xp, 9) +b{u- Up) , x^ (0) = 

Vs = (?Xs. (35) 
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According to (1161) . we have 

Xnew = Xp + Xs and y = yp + ys. (36) 

The strategy here is to assign the tracking (including rejection) subtask to the primary 
system (|34l) and the input-realization subtask to the secondary system (l35l) . It is clear from 
(I34l)-(l36l) that if the controller Up drives — )■ r in (l34l) and u drives — in (|35]) . 
then y — 7- r as t — 7- oo. The benefit brought by the additive state decomposition is that the 
controller u will not affect the tracking and rejection performance since the primary system 
(|34|) is independent of the secondary system (|35l) . Since the states Xp and Xs are unknown 
except for addition of them, namely Xmw, an observer is proposed to estimate Xp and Xg- 

Remark 5. Although the proposed additive state decomposition gives clear how to decom- 
pose a system, it still leaves a freedom to choose the primary system. By the additive state 
decomposition, the transformed system (l26l) can be also decomposed into a primary system 

Xp = Axp + bup, Xp (0) = 

Vp C Xp -\- dnew (37) 

and the derived secondary system 

/ ^t, Xp + Xs, — Axp + h {u^ — Up) , Xs (0) = Xq 



X c 



ys = c Xs (38) 

where A E ]R"X" is an arbitrary constant matrix. Therefore, there is an infinite number of 
decompositions. The principle here is to derive the secondary system with an equilibrium 
point close to zero as far as possible. If so, the problem for the secondary system is only 
a stabilization problem, which is easier compared with a tracking problem. In (|35]) . Xs = 
is an equilibrium point of Xs = f{t, Xp + Xs, 0) — f(t, Xp, 9), whereas in (1381) . = is not 
an equilibrium point of Xs = f {t,Xp + Xs,0) — Axp. This is why we choose the primary 
system as (|34|) not (l37l) . From the mention above, a good additive state decomposition often 
depends on a concrete problem. 

Theorem 3. Under Assumption 1, suppose that an observer is designed to estimate state 
Xp and Xs in (I34l)-(l35l) as follows: 

Xp = f{t, Xp, 0) + hup, Xp (0) = 

^new -Op. 
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Then Xp = Xp and Xs = Xs- 

Proof. Similar to the proof of Theorem 2. □ 

So far, we have transformed the original system to an uncertainty-free system, in which 
the new state and the new disturbance can be estimated directly. And then, decompose the 
transformed system into two independent subsystems in charge of corresponding subtasks. 
In the following, we are going to investigate the controller design with respect to the two 
decomposed subtasks respectively. 

B. Problem for Tracking ( including Rejection ) Subtask 
Problem 1. For design a controller 

Up U (t, Xp, r, C^jieio) (40) 

such that Vp ^ r + B {6r) ^ as t ^ oo, meanwhile keeping the state Xp bounded, where 

6r e R. 

Remark 6 (on Problem 1). Since Up = c^Xp + dnew, Problem 1 can be also considered 
to design Up such that c^Xp — )■ r — dnew Here, the difference between r and dnew should 
be clarified. The reference r is often known a priori, i.e., r (t + T) is known at the time t, 
where T > 0. Moreover, its derivative is often given or can be obtained by analytic methods. 
Whereas, the new disturbance dnew only can be obtained at the time t whose derivative only 
can be obtained by numerical methods. By recalling (1271) . the new disturbance dnew depends 
on the disturbance d, the parameters 9 and 6, the effect of unmodeled high-frequency gain 
namely ^, the state Xnew, and initial value xq. One way of reducing the complexity is to design 
an observer to estimate 9, and makes ^— T-^ast— T-oo.Asa result, the new disturbance dnew 
finally depends on d and ^ as t — t- cxo. In practice, low frequency band is often dominant in 
the reference signal and disturbance. Therefore, from a practical point of view, we can also 
modify the tracking target, namely r — dnew- For example, let r — dnew pass a low -pass filter 
to obtain its major component. If the major component of r — dnew belongs to a fixed family 
of functions of time, Problem 1 can also be considered as an output regulation problem [1211 . 

^B{5) = {xeR\ \x\<5}- the notation x{t) B (5) means min \\x (t) - y\\ 0; B (Si) + B [52) = 

ySl3(S) 

{x + y\xeB{Si),y eB{52)} 
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C. Problem for Input-Realization Subtask 

As shown in Fig.l, the input realization subtask aims to make ys — t- 0. Let us investigate 
the secondary system (l35l) . By Lemma 1, we have 

\\xsm</3{\\xs{to)ht-to) + ^\\b\\ sup \\u{s)-up{s)\\. (41) 

to<s<t 

This implies that ys ^ ^ \\h\\ \\c\\ B {5s) as u —i- Up + B (Ss), where Sg E M. It is noticed that u 
only can be realized by (|2TI) . Therefore, problem for input-realization subtask can be stated 
as follows: 

Problem 2. Given a signal Up, design a controller v = (t,Up) for ^271) such that 
u ^ Up + B {6s) as t ^ oo. 

This is also a tracking problem but with a saturation constraint. Here we give a solution 
to the Problem 2. The main difficult is how to handle the saturation in ((2T]) . Here, additive 
state decomposition will be used again. Taking (|2T]) as the original system, we choose the 
primary system as follows 

Zp = A^Zp + h^v 
Uzp = c^Zp + dzV (42) 

Then the secondary system is determined by the original system ((2TI) and the primary system 
(l42l) with the rule (fTSi) . and we can obtain that 

Zs = A^Zs + {cTa (v) - V) 

u^s = c^Zp + 4 (aa (v) - v) . (43) 

According to (fT6l) . we have z = Zp + Zg and u = Uzp + u^s- The benefit brought by the additive 
state decomposition is that the controller saturation will not affect the primary system (l42l) . 
Moreover, the controller v can be designed only based on the primary system (l42l) . where the 
controller v uses the state Zp not z. So, the strategy here is to design v = (t, Up) in (l42l) to 
drive u^p Up as t oo and neglect the secondary system (l43l) . Since v'^ {t, Up) is bounded, 
the state of the secondary system (l43l) will be bounded as well. If (Ja {v^ {t, Up)) —v^ {t, Up) — )■ 
as t — oo, then -> as t — )■ oo. Consequently, m — as t — )■ oo. For (l42l) . the 
transfer function from v to u^p is u^p (s) = C {s)v (s) . If C (s) is designed to be minimum 
phase, an easy way is to design v to be 

v{s) = C{s)up{s). (44) 
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The design will make the signal aa {v) close to the idea one, meanwhile maintaining the 
signal cTq (v) smooth as far as possible. By recalling (fTSl) . it will make the effect of the 
unmodeled high-frequency gain and the time delay ^ smaller. 

D. Controller Integration 

With the solutions of the two problems in hand, we can state 

Theorem 4. Under Assumptions 1-2, suppose i) Problems 1-2 are solved; ii) the controller 
for system ([T]) (or (|26l) ) is designed as 
Observer: 

&new f(j^i ^newi ^) ~l~ hu^ Xnew (0) 0, 

Xp = f{t, Xp, 9) + hUp, Xp (0) = 0, 

dnew y C Xnew (45) 

Controller: 

Up U (t, Xpj r, dnew) i ^ ^ (tj ^p) 

z = A^z + h^aa (v) ,u = c^z + d^cXa (v) (46) 

Then the output of system ([T]) (or (|26l) ) satisfies that y r + B {5r + 'j \\b\\ \\c\\ 5s) as t — )■ oo, 
meanwhile keeping all states bounded. In particular, if 5r + 5s = 0, then the output in system 
© (or (|26l) ) satisfies that y — )■ r as t — )■ cxd. 

Proof. It is easy to follow the proof in Theorems 2-3 that the observer (l45l) will make 



i"new — •^newt dnew — dnewi — -^p' (4V) 

Suppose that Problem 1 is solved. By (l40l) and (l47l) . the controller Up = u'^{t, Xp, r, dnew) can 
drive yp r + B {5r) as t — )■ oo in (|34|) . Suppose that Problem 2 is solved. By (l47l) . the 
controller v = (t, Wp) can drive u Up + B {5s) as t — oo in (l35l) . Further by (HTI) . we 
have ys^ Bi-f \\b\\ \\c\\ 5s) . Since y = yp + ys, we have y -> r + i3 (5^ + 7 ||c|| 5,) . □ 

Example 5 {Rohrs' Example, Example 3 Continued). According to (|34l) . the primary system 
of linear system (l32l) can be rewritten as follows: 

Xp = -{3 + 9)xp + 2up 

y -^p ~l~ dn^yj 
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Design Up = |[(2 + 9)xp + r + r — dnew — dnew]- Then the system above becomes Cp = —Cp, 
where Cp = i/p — r. Therefore, ?/p — t- r as t — )• oo. According to (|44l) . v is designed as 

(t, Up) = 2up + Up. Here iip and (i„e«> are approximated by dnew ~ '^""'^(oif+T'^'^e'" (^)) ^^'^ 
iip ~ (s)), respectively. Suppose ^ = and given r = 0.5 and r = 0.5sin(0.2t) , 

respectively. Driven by the resulting controller (l46l) . the simulation result is shown in Fig. 2. 

Example 6 (Nonlinear, Example 4 Continued). According to (|34l) . the primary system of 
nonlinear system (|33l) can be rewritten as follows: 

Xp = -Xp - (1 + + Up, Xp (0) = 

Up Xp -\- dnew 

Design Up = (1 + + r + r — (i„e«) — dnew Then the system above becomes Cp = 
—Cp, where Cp = Up — r. Therefore, — )■ r as t — > oo. According to (l44l) . f'' (t, Up) is 
designed as v'^ (t, Up) = 2up + Up. Here the derivative of Up and dnew are approximated by 

dnew ~ '^"^(aTfTT'^"^"' ^'^^ ~ '^"^(alf+T^p^'^))' respectively. Suppose ^ = and 
given r = 0.5 and r = 0.5sin(0.2t) , respectively. Driven by the resulting controller (l46l) . the 
simulation result is shown in Fig. 3. 

Remark 7. The derivative of dnew and Up can be also obtained by differentiator technique 
[|23]| . [|24l . It is interesting to note that 9 is different from 9, but y ^ r with a good tracking 
accuracy. This is one major advantage of this proposed control scheme. Moreover, all the 
unknown parts such as 9, d and the effect of H (s) e"*"" are treated as a lumped disturbance 
dnew This can explain why the proposed scheme can handle many uncertainties. 

V. Two-Cart Example 

The two-cart mass-spring-damper example was originally proposed as a benchmark prob- 
lem for robust control design [|T2l . [l22l . Next, we will revisit the two-cart example by the 
proposed control scheme. 

The two-cart system is shown in Fig. 4. The states xi(t) and X2{t) represent the absolute 
position of the two carts, whose masses are mi and m2 respectively; ^1,^2 are the spring 
constants, and bi,b2 are the damping coefficients; d(t) is a disturbance force acting on the 
mass m2; u{t) is the control force subject to an unmodeled high-frequency gain and a time 
delay, which acts upon the mass mi. The parameter mi = 1 is known, whereas the following 
parameters m2 = 2, A;i = 0.8, k2 = 0.5, bi = 1.3, 62 = 0.9 are assumed unknown. The un- 
modeled high-frequency gain and a time delay is assumed to be H (s) e~'^^ = ^2^3o^g%29 ^~°^^- 
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The disturbance force ({t) is modeled as a first-order (colored) stochastic process generated 
by driving a low-pass filter with continuous-time white noise s{s), with zero-mean and unit 



intensity, i.e. S = 1, as follows ( (s) 



0.1 
s+O.l 



eis 



The overall state-space representation is formulated into ([T]) as follows: 
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The objective here is to design a tracking controller u such that y r with a good tracking 
accuracy. Since the dampers will always consume the energy untile the two carts are at rest, 
it can be concluded that the two-cart system (a physical system) is stable for any 6 in the real 
world. This implies that, for any given 9 E Vl, there exist positive definite matrices P G M"^" 
and Q G M"''" such that PA (9) + A^ (9) P < -Q, \Jx G M", where 9. represents the set that 
any 9 in the real world. Assumption 1 is satisfied. Let C (s) = Then Assumption 2 is 
satisfied with eh = 0.12 and = 1. 

According to input redefinition above, the two-cart system (l48l) can be rewritten as follows: 

x = A{9)x + bu + {d + ^) 

T 

y = c x 

where sup^>o 1^ (^)l ^ 0.17a and u is generated by 

z = -0.5z + 0.5(Ta (v) , 



u 
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Then, according to (|26l) . the uncertainty-free transformed system of (l48l) is 



?/ C Xfiew ~l~ dnew (49) 

Where d„e^ = c'^e^(^)*xo + / c'^e^('^^'-'^b[d [s) + ^s) + A [6) x^eu, {s) - A{e)xneu, {s)]ds. 

Jo 

According to (l34l) . the primary system of (|49l ) can be rewritten as follows: 

Xp = A{9)xp + bup, Xp (0) = 

Up C Xp -\- dnew (50) 

The transfer function from Up to yp in (l50l) is G^^^ (s) , which is a minimum phase. Thus, Up 
can be designed as Up (s) = Gy^ (s) (r — dnew) (s) , which can drive yp — )■ r. The Problem 
1 is sloved. Furthermore, according to (l44l) . redefined input v is designed as t>'' {t, Up) = 
C^^ (C^^ (s) Cy^^ (s) (r — dnew) (s)) . To realize the control, {t, Up) is approximated to be 

{t, Up) = {Q (s) G-' (s) Gy^ (s) (r - d^ew) (s)) . (51) 

where Q (s) is a fifth-order low-pass filter to make the compensator physically realizable (the 
order of denominator is greater than or equal to that of numerator). For simplicity, Q (s) is 

chosen to be Q is) = / \ r here. The Problem 2 is sloved. Therefore, according to 

(l45])-(l46l). the controller for the two-cart system is designed as follows: 

^new A(^0)Xfiew ~t~ bu^ -^new (0) d^ew V ^ ^new 

z = -0.5z + 0.5aa{v'' {t,Up)) ,u = z (52) 

where {t, Up) is given by (|5TI) and a is chosen to be 1 here. 

To shown the effectiveness, the proposed controller (l52l) is applied to three cases: 

Case 1:9 = 6 

Case 2: ^ = [ 1 l l 0.9 1.5 1 Y 

Case 3:e = e=[l l l 0.9 1.5 1 V- 
Case 1 implies the parameters are known exactly. Case 2 implies the parameters are unknown. 
While, Case 3 implies the parameters are changed to be a specified one. The simulations are 
shown in Figs. 5-7. The proposed controller achieves a good tracking accuracy. Moreover, it 
is seen that the response in Cases 2-3 is faster than that in Case 1 . And, the tracking accuracy 
in Cases 1,3 is better than that in Case 2. So, Case 2 is a tradeoff between Case 1 and Case 
3. 



January 26, 2012 



DRAFT 



20 



Remark 8. The simulations show that the proposed controller can handle the case that the 
estimate parameters are different from the true parameters. Moreover, the response is similar 
to that of the model with the estimate parameters. This implies that the proposed controller in 
fact achieves the results similar to model reference adaptive control. However, unlike model 
reference adaptive control, unknown parameters are not estimated and changed directly. 

Remark 9. If the considered system is parameterized but there exist many uncertain 
parameters, then an adaptive control often needs corresponding number of estimators, i.e., 
corresponding number of integrators. This will cause parameters converging to true values 
with very slow rate or cannot converge to true values if without persistent excitation. Whereas, 
in the proposed control, five uncertain parameters and disturbance are lumped into the 
disturbance dnew, which can be estimated directly. 

VI. Conclusions 

Output tracking control for a class of uncertain systems subject to unmodeled dynamics 
and time delay at input is considered. Our main contribution lies on the presentation of a new 
decomposition scheme, named additive state decomposition, which not only transforms the 
uncertain system to be an uncertainty-free system but also simplifies the controller design. The 
proposed control scheme is with the following two salient features, (i) The proposed control 
scheme can handle mismatching uncertainties and mismatching disturbance. Moreover, it 
can achieve a good tracking performance without exact parameters, (ii) The proposed control 
scheme has considered many uncertainties. In the presence of these uncertainties, the closed- 
loop system is still stable when incorporating the proposed controller. Three benchmark 
examples are given to show the effectiveness of the proposed control scheme. 
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Fig. 2. Output of the Rohrs' example system 
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Fig. 3. Output of the nonhnear system 
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Fig. 4. The two-cart system 
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Fig. 5. Output of the two-cart system in Case 1 
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Fig. 6. Output of the two-cart system in Case 2 
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Fig. 7. Output of the two-cart system in Case 3 
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